XMM detection of the supernova remnant G337.2+0.1 
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\D '. ABSTRACT 
O 

Oh 1 We report the first XMM detection of the SNR candidate G337.2+0.1 

(=AX J1635.9— 4719). The object shows centrally filled and diffuse X-ray emis- 
sion. The emission peaks in the hard 3.0—10.0 keV band. A spatially resolved 
spectral study confirms that the column density of the central part of the SNR is 
about Nh ~5.9(±1.5) x 10 22 cm -2 and its X-ray spectrum is well represented by 
^ a single power-law with a photon index T=0.96±0.56. The non-detection of line 

emission in the central spectrum is consistent with synchrotron radiation from 
a population of relativistic electrons. Detailed spectral analysis indicates that 
the outer region is highly absorbed and quite softer than the inner region, with 
Ntf ~16.2(±5.2)xl0 22 cm" 2 and kT=4.4(±2.8) keV. Such characteristics are al- 
ready observed in other X-ray plerions. Based on the morphological and spectral 
X-ray information, we confirm the SNR nature of G337.2+0.1, and suggest that 
the central region of the source is a pulsar wind nebula (PWN), originated by an 
energetic though yet undetected pulsar, that is currently losing 
of ~ 10 36 erg s -1 . 



Subject headings: X-rays: individual (AX J1635.9— 4719) — (ISM:) supernova 
remnants — ISM: individual (G337.2+0.1) — X-rays: ISM — radiation mecha- 
nisms: non-thermal 
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1. Introduction 

Supernova remnants (SNRs) of the Crab-like or plerionic class are objects characterized 
in the radio band by a compact, filled-center morphology with a relatively flat spectral index 
(Slane et al. 2000). In the X-ray band they present non-thermal spectra, characteristic of 
synchrotron emission. The non-thermal spectrum, in some cases, can even reach the very 
high-energy gamma-ray region, like in MSH 15-52, G18.0-07, and Vela X (see Aharonian et 
al. 2005, 2006a). It is widely believed that plerions are powered by the loss of rotational 
energy from energetic pulsars, although clear evidence of the presence of these objects is 
often lacking. The pulsar wind forms a nebula inside the SNR (the PWN), where relativistic 
particles can be efficiently accelerated producing synchrotron radiation that yields the typical 
morphologies observed at radio and X-rays. 

The new generation of X-ray instruments with improved sensitivity and spatial/spectral 
resolution like ASCA, XMM, CHANDRA and INTEGRAL allows to detect and study not 
only distant energetic objects with these characteristics, but also those located in regions of 
high density. 

Recently, Combi et al. (2005) have presented evidence supporting a SNR origin for the 
radio source G337.2+0.1. The object was discovered by the MOST galactic plane survey at 
843-MHz towards the Norma spiral arm (Whiteoak & Green 1996), and later, detected with 
the ASCA telescope during a survey of part of the galactic plane (Sugizaki et al. 2001). Its 
integrated flux was reported to be ~ 1.2xl0 -12 erg cm -2 s" 1 in the 0.7-10 keV energy band. 
The photon index is poorly constrained: r=2.8^ 2 ;g, and the absorption column density of the 
best-fit model yielded N^=15±g 5 x 10 22 cm" 2 (Combi et al. 2005). A thorough study of its 
radio (continuum and line) and X-ray properties shows that the emission from the source is 
consistent with what is expected for a young SNR located at a distance d ~14 kpc (Combi et 
al. 2005). More recently, G337.2+0.1 has been suggested as the potential counterpart of the 
high-energy gamma-ray source HESS J1634-472 (Aharonian et al. 2006b): this possibility 
needs further confirmation. Throughout this paper, we adopt 14 kpc as the distance to 
G337. 2+0.1 (hence, 1' corresponds to 4 pc) . 

In this Letter, we present the first XMM-Newton observations of the SNR candidate 
G337. 2+0.1. Based on its X-ray properties we are able to confirm that this object is a non- 
thermal SNR with hard, featureless power-law spectrum and, possibly, a PWN originated by 
a non-detected energetic pulsar. In Section 2, we present the X-ray observations and data 
reduction. X-ray analysis and results are presented in Section 3, and in Section 4 we provide 
a brief discussion and a summary. 
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2. X-ray Observations and data reduction 

The SNR candidate G337.2+0.1 has been marginally observed on February 2004 by 
the XMM-NewtonX.-ray satellite in two separated observations (Obs-Id. 0204500201 and 
0204500301). Both observations were centered towards the source IGR J16358-4726 (aj2ooo.o= 
16 h 35 m 53 ! !820, 5 J20 oo.o=-47 25'41 / . / 10), and were acquired with the EPIC MOS (Turner et al. 
2001) and EPIC PN (Striider et al. 2001) cameras. Observations were taken with a n thin n 
filter, and in the full frame (FF) imaging mode. Temporal resolution is 2.5 s and 200 ms for 
the MOS and PN CCDs, respectively 1 . Observations were obtained from the XMM- Newton 
Science Archive (XSA) 2 , and raw EPIC data were calibrated using the last version of the 
Standard Analysis System (SAS) 3 . To create images, spectra and light curves, we selected 
events with flag 0, and patterns< 12 and 4 for MOS and PN cameras, respectively. 

The off-angle of G337. 2+0.1 respect to the center of the observation is about ~6.67 
arcmin, which implies a reduction of the XMM-Newton effective area of about 12%. The 
net exposure times of the observations are 34.7 (Obs-Id. 0204500201) and 32.6 (Obs-Id. 
0204500301) ks. Unfortunately, the first observation was affected by a high and variable 
soft proton background level (Lumb et al. 2002), whereas the second one (Id. 0204500301) 
is unaffected by background fluctuations. We derive Good Time Intervals (GTI) by the 
accumulation of background light-curves in the 10-15 keV energy band, which leads to a 
reduction of ~87% in the net exposure time of the Obs-Id. 0204500201. In order to avoid 
contamination for high background patterns hereafter our analysis concerns only to the 
observation 0204500301. The number of detected counts in the 0.5-2.5 and 2.5-10.0 keV 
energy bands are 117/121/315 and 431/403/1154 for the MOS1, MOS2 and PN cameras, 
respectively. Finally, at the SNR G337.2+0.1 EPIC-PN position there is a CCD gap in the 
X-ray image, leading us to ignore these data only for the X-ray image analysis section, but 
they are included for the rest of our study. 



1 http://xmm.vilspa.esa.es/external/xmm_user_support/documentation/ 
2 http://xmm. vilspa.esa.es/xsa/ 

3 http://xmm. vilspa.esa.es/external/xmm_sw_cal/sas.shtml 
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3. X-ray analysis of G337.2+0.1 
3.1. Image 

The coordinates of the SNR G337. 2+0.1 were defined at the position where X-ray emis- 
sion peaks (a J20 oo.o = 16 h 35 m 54!95, 5 J20 oo.o = -47°19'02"2), being the errors in R.A. and 
DEC of e a = ±251 and e$ = ±3'.'3, respectively (at the 90% of confidence). This position 
agrees well with the previous estimate of the radio position (Combi et al. 2005), but differs 
in ~ 50 arcsec from the ASCA coordinates. Because of the poor spatial resolution (~2.9 
arcmin) the ASCA telescope alone is not conclusive in resolving the SNR G337. 2+0.1. For- 
tunately, the availability of the 0.5-10.0 keV X-ray observations from XMM-Newton , the 
largest X-ray telescope so far (1480 cm 2 @1.5 keV), allows to investigate with success the 
X-ray nature of this source. In fact, XMM-Newton has ~40 times more spatial resolution 
than ASCA, whereas the sensitivity limit is at least ~3 times better. 

Figure 1 shows the X-ray image of the SNR G337. 2+0.1 in the 0.5-10.0 keV energy 
band. The image does not reveal a typical rim-brightened outer SNR shell, so the overall 
size of the diffuse X-ray emission is uncertain. 

We use the clean event files to generate MOS1 ans MOS2 images in the energy band 
[0.5-10] keV with a spatial binning of 4.35 arcsec per pixel. In order to increase the signal- 
to- noise (S/N) ratio, we use the emosaic SAS task to merge together the two images. The 
corresponding set of exposure maps for each camera has been prepared to account for spa- 
tial quantum efficiency and mirror vignetting by running the SAS task eexmap. Exposure 
vignetting corrections were performed by dividing the superposed count image by the corre- 
sponding superposed exposure maps. We adaptively smoothed this image to a S/N ratio of 
10 using the SAS task asmooth (see Figure 1). Plotted contours correspond to 1, 2, 3 and 4 
sigma levels over the mean background flux of the image (~6xl0~ 5 ph/cm 2 /s). 

Finally, we are able to investigate the spatial extent of G337.2+0.1. Initially, this was 
performed by Combi et al. (2005) using ASCA data. However, the low resolution of this 
image led these authors to consider a radial analysis at large angular distance (~ 10 arcmin). 
Thanks to the high resolution XMM-Newton image, we found at least 19 weak point-like 
X-ray sources inside this radius, making the ASCA studies of G337.2+0.1 to be hardly biased 
by source contamination effects. According to the image presented in Figure 1, G337.2+0.1 
does not extend farther than 1.5 arcmin from the central peak. We also compare G337.2+0.1 
spatial extent with that produced by a point source placed at a similar off-axis (~ 6 arcmin) 
position. In Figure 2 we show that the SNR G337. 2+0.1 has an extension ~3.5 times larger 
than what is expected for a point-like source. 
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3.2. Spectral analysis 

For the spectral analysis we used MOS and PN data. It was performed using the xspec 
package (Arnaud 1996). Since the statistics of the source is not complete enough to perform a 
spatial-spectral analysis, we extracted X-ray photon events from only three different regions: 
i) a circular region of 50 arcsec; ii) a circular region of only 12 arcsec that accounts for the 
central source observed in the image; Hi) an annulus for the extended emission of inner-outer 
radii of 12-50 arcsec. The background region was taken from a nearby blank region in the 
neighborhood of the source. Ancillary response files (ARFs) and redistribution matrix files 
(RMFs) were calculated. All spectra are grouped with a minimum of 16 counts per bin. 
The background-subtracted spectra of the MOS and PN data (the upper line) are shown in 
Figure 3. We checked for possible background contamination in our spectra, inconsistency 
in the extraction of local backgrounds and differences in background corrected spectra of 
SNR G337. 2+0.1. We note excesses in spectra, essentially at energies between 7.5 to 8 keV. 
However, such features have a low statistical significance (~ 1 to 1.5 sigma) and are related 
to fluorescence lines in the background spectrum of the XMM-Newton (e.g. De Luca & 
Molendi 2004). 

Our analysis of the XMM-Newton EPIC spectra was essentially performed using a 
single non-thermal model, described by a simple Power-Law emission model. We also fit the 
spectra by a thermal emission model (apec, Brickhouse 2003). Both models were affected 
by an absorption ISM component (WABS, Morrison & McCammon 1983). The goodness of 
the model fit was derived according to the x 2 -test statistics. The best fitting parameters of 
the models are shown in Table 1, and the errors quoted are 90% confidence limits. 

According to the results presented in Table 1, the central part of the SNR appears quite 
harder (r ~ 0.96) than the outer one (r ~ 2.38). We suggest that the most reasonable 
interpretation of observed emission from the central part of SNR is synchrotron radiation 
from relativistic electrons accelerated in the vicinity of the central source of the SNR. 

The softening of the spectrum toward the outer regions of the nebula is a well known 
effect which has been seen in other X-ray plerions (e.g., GO. 9+0.1, Porquet et al. 2003; 3C58, 
Torii et al. 2000; G21. 5-0.9, Slane et al. 2000). 

To get a statistical assessment of the X-ray variability of the SNR G337.2+0.1, we use the 
32.6-ksec EPIC-PN observation to compare the time arrival distribution of source photons 
by means of the Kolmogorov-Smirnov (KS) test (Press et al. 1992). We use an extraction 
region centered in the SNR with a radii of 50 arcsec. The total number of photons is 1470. 
We see no significant pulsed signal with a period greater than twice the read-out time of the 
EPIC-PN camera in the FF mode, which corresponds to a Nyquist limit of 400 ms. 
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4. Discussion 

The X-ray morphology of G337.2+0.1 shows a centrally peaked emission, surrounded by 
a diffuse X-ray nebula. The lightcurve of the object does not show any significant flux vari- 
ability above 0.4 seg, implying that at first glance, a pulsar origin for the central contribution 
could be ruled out. However, a detailed spectral analysis indicates that the outer region is 
softer than the inner region, a phenomenon observed previously in several X-ray plerions 
with PWN (e.g., IC443, Bocchino & Bykov 2001). A spectral analysis of the central com- 
ponent of the SNR shows that the X-ray spectrum is well represented by a single power-law 
with a photon index T=0.96±0.56, a value similar to that of objects powered by an energetic 
pulsar (Gotthelf 2003). Moreover, the non-detection of line emission in this spectrum is 
consistent with synchrotron radiation from a population of relativistic electrons. These facts 
suggest a non-thermal origin for the X-ray emission. We therefore conclude that the system 
G337. 2+0.1/ AX J1635.9— 4719 is a non-thermal SNR with, possibly, a non-detected pulsar. 

Possible reasons for the non-detection of a pulsar inside the SNR are a short rotation 
period (less than 400 ms) or unfavorable geometrical conditions. The presence of a pulsar is 
suggested by the central X-ray peak found inside G337.2+0.1. In what follows we explore the 
possibility that there exists a hidden pulsar-powered component (plerion) within the SNR. 
Using the empirical formula derived by Seward & Wang (1988), log L x (erg s _1 ) =1.39 log E — 
16.6, where L x is the X-ray luminosity of the plerion in the 0.2-4 keV band, we can make an 
estimate of the spin-down luminosity of the pulsar (see, also, Becker & Triimper 1997). Using 
the X-ray flux of the compact source and its nebula, Fx(0.2-4 keV)=4.9±1.7x 10~ 13 erg s _1 
cm -2 , we get Lx=l.lxl0 34 erg s -1 (unabsorbed). This implies a spin-down luminosity of 
E ~2.5xl0 36 erg s -1 , and a period of P >0.08 (£3-^33) 1//2 (s), where E^ is the spin-down 
luminosity in units of 10 38 erg s _1 , and t 3 is the pulsar age in units of 10 3 years. In order 
to compare this result with those obtained with other empirical relations between the X-ray 
luminosity and the rate of the spin-down energy loss, we have used the Becker & Triimper 
(1997) and Possenti et al. (2002) equations. In the first case, (taking into account only the X- 
ray flux of the point source in the 0.1-2.4 keV band) the spin-down luminosity is E ~3xl0 35 
erg s _1 , a factor 9 lower than the value obtained with the Seward & Wang (1988) relation. 
In the second case (using the X-ray flux of the compact source and its nebula in the 2-10 keV 
band), we got E ~8xl0 36 erg s -1 , a value that is a factor 3 higher than the value obtained 
with the Seward & Wang (1988) relation. On average, a value ~ 10 36 erg s" 1 seems to be 
reasonable. If we assume a pulsar period of less than 0.4 s, we obtain an upper limit for the 
age of the pulsar of t < 1000 years. 

We have seen that G337. 2+0.1 does not show a rim-brightened outer SNR shell. This 
could be the result of the absorption of the soft thermal emission from the forward shock 
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by the very high absorbing column density. Other sources like Crab, G21.5-0.9 (Slane et al. 
2000) and 3C58 (Torii et al. 2000) have weak or absent X-ray rims and all are powered by 
young X-ray pulsars (Murray et al. 2002; Camilo et al. 2006). 

It could be interesting to compare the characteristics of G337. 2+0.1 with 3C58. The 
X-ray luminosities, between 0.5 and 10.0 keV, are ~ 4.8 x 10 34 erg s^ 1 and ~ 2.4 x 10 34 erg 
s -1 , respectively. The radio luminosities, at 1 GHz, are ~ 3 x 10 32 erg s" 1 and ~4x 10 32 
erg s -1 . We see, then, that both sources are quite similar. We notice that the estimated age 
of 3C58 is ~ 800 yr, a value based on the association of the SNR with the supernova 1181 
(Stephenson & Green 2002). The nature of this SNR is discussed by Camilo et al. (2006), 
who, however, argue against the association with this historical SN (Bietenholz et al. 2001). 
The most significant difference seems to be the absence of any thermal component in the 
case of G337. 2+0.1. If we compare with the Crab, on the contrary, we see that the Crab 
pulsar is injecting around two orders of magnitude more energy per time unit in the nebula 
than G337. 2+0.1. The spin-down luminosity inferred for the pulsar from the new X-ray 
data sets an upper limit to the energy available for high-energy cooling channels like inverse 
Compton scattering and proton-proton interactions. The luminosity of the nearby HESS 
source J1634-472 (E > 1 TeV), if it is located at the same distance inferred for G337.2+0.1, 
would be L 7 ~ 7 x 10 34 erg s -1 . So, a physical association would be possible only if ~ 7 % 
of the spin-down luminosity is converted in high-energy 7-rays. 

Complementary studies of the PWN scenario will involve high-resolution X-ray obser- 
vations with Chandra satellite, and radio observations with ATCA, to allow the comparison 
of the X-ray spectrum and morphology with those at the radio band. GLAST observations 
could reveal a GeV 7-ray source if the proposed association with HESS J1634-472 is correct. 
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Table 1: X-ray spectral parameters of the SNR G337.2+0.1 



Region: 


All 


Inner 


Outer 


Model: 


Power-law 


APEC 


Power-law 


APEC(f) 


Power-Law 


APEC 


(cm" 2 ) in [10 22 ] 


11.49±2.73 


11.17±2.18 


5.9±1.52 


6.45±2.31 


16.21±5.62 


16.26±5.25 


T | kT(keV) 


1.82±0.45 


10.01±6.49 


0.96±0.56 


(64±64) 


2.38±0.78 


4.38±2.79 


Abundance* 




<0.05 




<0.05 




<0.05 


Norm, in [xlO -4 ] 


4.91±1.93 


8.28±2.83 


0.55±0.51 


(0.83±2.4) 


4.01±1.44 


6.7±2.81 


Flux (ergs- 1 cm _2 ) 


ll-7 8 2 8 9 


9.81& 1 


1 3 21 


i il.4 
i - i 0.9 


o o9.9 
°- Z 7.2 


7 49.3 


Reduced xl 


1.13 


1.14 


1.10 


1.25 


1.09 


1.06 



Notes: Flux is absorption-corrected in units of 10 13 (erg s 1 cm 2 ) calculated 
in the (0.5-10.0 keV) energy range. Normalization was calculated according to: 
10 14 /47rr 2 (l — z) 2 j n e n H dV. f APEC thermal model could not fit the spectra 
in a consistent way, yielding very ill-constrained parameters as solution. The 
non-thermal (power-law) fit seems to be the most representative emission model 
for the central part of the SNR. f The abundance was adopted from Anders &c 
Grevesse (1989) and left as free parameter along all our fits, but its value is 
hardly affected by the low statistic of the spectra. 
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Fig. 1. — X-ray image of the G337. 2+0.1 in the 0.5-10.0 keV band. Typical local back- 
ground level of the smoothed image is ~0.94(±0.17) counts px _1 . Contours show the level 
of 1.6 (~1<t), 3.2 (~2<r), 4.7 (~3cr) and 6.3 (~4<r) photons px -1 , from outer to inner curves, 
respectively. 
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Fig. 2. — Radial profile of the smoothed X-ray image of SNR G337. 2+0.1 (filled circles). 
Open squares refer to the radial profile of an observed point source located at roughly the 
same off-axis position. 
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Fig. 3.— X-ray spectra of the SNR G337.2+0.1 in the 0.5-10.0 keV band. In all figures PN 
data are the upper line. Left: Spectrum extracted on the whole SNR. Center: Spectrum of 
the compact source observed in the center of SNR G337.2+0.1. Right: Outer spectrum that 
excludes the central contribution. 



